a=1,b=100

xa=-2:0.1:2; ya=-1:0.1:3;

[x,y]=meshgrid(xa,ya);

z=(a-x).^2+b*(y-x.^2).^2;

mesh(x,y,z)

xlabel('x');ylabel('y');zlabel('z');